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Abstract: The paper studies absolute stability of neutral differential nonlinear systems 
x(t) = Ax (t) + Bx (t - T) + Dx (t — 1) + bf(o(t), o(t) = cx(t), t 20 


where x is an unknown vector, A, B and D are constant matrices, b and c are column constant vectors, T > 0 
is a constant delay and f is a Lurie-type nonlinear function satisfying Lipschitz condition. Absolute stability 
is analyzed by a general Lyapunov-Krasovskii functional with the results compared with those previously 
known. 
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1 Introduction 


The paper considers absolute stability of a system of nonlinear differential equations of neutral type 
x(t) = Ax(t) + Bx(t - T) + DX(t- T) + bf(o(), t20 (1.1) 


where x = (x1,..., Xn)”: [-7, œ) > R”, A, Band D are nxn constant matrices, b is an n-dimensional column 
constant vector, T > 0 is a constant delay, f: IR > R, f (0) = 0, is a Lurie-type nonlinear function satisfying 
Lipschitz condition, 

olt) :- c'x(t) (1.2) 


and c is a column n-dimensional constant vector where the superscript T denotes the transpose of a vector 
(or a matrix). 


1.1 On absolute stability 


The problem of absolute stability was formulated in 1944 by Lurie and Postnikov in [39]. In [1, 8], general- 
izations of the problem are carried out, and the well-known Aizermann hypothesis [1, 49] formulated (some 
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scholars still adhere to referring to the problems associated with absolute stability as Lurie nonlinear control 
system ones). The initial period of studying the theory of absolute stability is associated with the name of its 
founder [39, 40], who tried to solve the problem with the direct Lyapunov method by means of a Lyapunov 
function of the “quadratic form plus the integral of nonlinearity” type. This classical method for investigat- 
ing absolute stability using Lyapunov's functions was used by many authors [1, 40-42, 60]. Since the mid 
1960's the intensive development in the theory of absolute stability has included the so-called frequency 
method, initiated by Popov [21, 22, 43, 45]. In addition to the direct use of so called Popov's criterion for 
the study of absolute stability, special methods for its application are sought and analogous frequency cri- 
teria are formulated [2, 25]. An in-depth analysis of both of the above methods, a study of the relationship 
between them, and a number of new original approaches are contained in the numerous publications by 
Yakubovich and his followers [21, 22, 34, 59]. When constructing mathematical models of real systems for a 
more precise description of the operation of systems, there is a need to consider equations with delay and, 
in general, neutral type equations (for a recent overview we refer to [48]). Unfortunately, time-delay phe- 
nomenon can frequently degrade performance of control systems or even can lead to instability. Therefore, 
the research of stability in general and, particularly, absolute stability of time-delay Lurie systems is of im- 
portant significance [10, 11, 17, 19, 26, 28-30, 32, 35, 47]. In [31], Krasovskii replaces Lyapunov functions by 
what is now called Lyapunov-Krasovskii functionals, which find ample applications in setting out the con- 
ditions of stability. Numerous ingenious constructions have been published of such functionals, often using 
the so-called Razumihin condition [50]. In particular, we refer to [5, 9, 20, 46, 51—53, 56, 58]. Surveys and 
overviews of the methods and results related to absolute stability achieved since the birth of the theory in 
1944 can be found, for example, in [2, 21, 32-34]). Note that a positivity-based approach, closely related to 
the concept of the Lyapunov functions, which uses the nonnegativity of entries of the Cauchy matrices, is 
developed in [3, 4, 6, 7, 12, 13, 15]. The above methods are widely used in investigating the stability of neutral 
differential equations. For recent results we also refer to [14, 16, 24, 36, 54, 55, 57] and to the references therein. 


1.2 Preliminaries, the problem considered, and the paper structure 


Throughout the paper we use the following restriction on the nonlinearity of f (not mentioning this explicitly 
when formulating results). Namely, assume that there exists a k > O such that 


kw - f(m)| f(w) » 0, Vw #0. (1.3) 
[ ] 


Such a nonlinearity is often said to be of a Lurie type or, based on a geometrical interpretation of (1.3), we 
also say that f satisfies a sectorial-condition. 
Let o: [-T, 0] > R” bea continuously differentiable vector-function. In addition to system (1.1), consider 
an initial condition 
x=), x(t- e(t. tel-r.ol. (1.4) 


A solution of the problem (1.1), (1.4) is understood within the meaning of the following definition. 


Definition 1. A function x: [-1T, œ) > R” is called a solution of (1.1), (1.4) if it is continuous on [-T, co), con- 
tinuously differentiable on [-T, œ) V (Ts, s = 0,1,...}, satisfying (1.1) on [0, œ) V (Ts, s = 0,1,...}, and 
satisfying (1.4). 


The existence of a solution of the problem (1.1), (1.4) on [-T, co) within the meaning of Definition 1, its unique- 
ness and continuous dependence on initial values (1.4) are the consequences of the linearity of the first three 
terms in the right-hand side of (1.1), the Lipschitzeanity of f and the sectorial condition (1.3). 

For a vector y = (y1, ..., yn)! : [-T, eo) > R”, the following norms are used in the paper 


n 1/2 
LGE (3-70) , if te[-r, e9), 
i=1 
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lIVller = max , {|y(s +O|}, if te[0,e9), 


i 1/2 


lees ( f es totas |. if ceto. 


eT 


where c > Oisaconstant. If B is a symmetric positive-definite matrix, denote by Amax (2), Amin(B) its maximal 
and minimal eigenvalues. We also use a norm of a square matrix A, induced by above Euclidean norm of a 
vector, computed by formula 


[Al = (AmasCAT4)) (1.5) 


Listed below are the definitions of stability used. 


Definition 2. The zero solution of system (1.1) is exponentially stable in the metric C? if, for the solution 
x: [-T, eo) > R” of (1.1), (1.4) we have 


Ix(£) < Fog, ge", t0, (1.6) 


where 3) is a positive definite functional and qo > O is a constant. 


Definition 3. The zero solution of system (1.1) is exponentially stable in the metric Ct, if it is stable in the metric 
C? and, for t € [0, œ) V(15,$ = 0, 1, ... 


Ix (5| « Flp, gem, (17) 


where 31 is a positive definite functional and q4 > O is a constant. 


Definition 4. The system (1.1) is said to be absolutely stable if its zero solution is exponentially stable in the 
metric C! for an arbitrary function f, satisfying (1.3). 


In the present paper we solve the problem of the absolute stability of system (1.1) within the meaning of 
Definition 4. The paper is structured as follows. First, the exponential C stability of (1.1) is studied in part 2 
with the main result of this part being Theorem 1. Its proof is carried out by the Lyapunov-Krasovskii method 
of functionals using a special functional V. Such an approach needs the positivity of a matrix (the matrix S 
below) appearing in the estimate of the derivative of the functional V along solutions of (1.1). Second, the 
exponential C! stability of (1.1) is studied in part 3. The proof of the main result of this part (Theorem 2) is 
based on the possibility to exclude the delayed derivative term Dx(t — T) from (1.1), considering then (in an 
admissible domain) only the delayed system derived instead (Lemma 1). Then, the final statement (Theorem 3) 
is formulated. The results are illustrated by an example in part 4. The paper is concluded by part 5 with 
comments and conclusions concerning the approach used. Relations with the investigations of other authors 
are discussed as well. In the formulations of the theorems on C? or C! exponential stability, functionals 3; 
and constants ri, i = O, 1 are defined explicitly with their constructions following from the methods of proofs. 


2 Exponential C? Stability 


Below, we give conditions for C? exponential stability of (1.1). From the proof, an estimate of the convergence 
of solutions follows as well. To solve the problem of absolute stability within the meaning of Definition 4 for 
system (1.1) under restriction (1.3), the Lyapunov-Krasovskii functional (2.1) will be used of the quadratic type 
with respect to both the current coordinates and their derivatives, 


t t o(t) 


V(x(£) = x (t) Hx(t) + I EdT (s) Gax(s)ds + I e $7957 (s) GX (s) ds +B I f()ds. 23) 
0 


tr t-T 
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In (2.1), t > 0, x isa solution of system (1.1), H, G1, G2 are positive-definite symmetric matrices, çı > 0, ¢2 > 0 
and f > 0. The condition (1.3) and the well-known properties of quadratic forms imply the two-sided estimate 
of functional (2.1) 


Amin (H) |X (| + Amin(G1)|IX|ler,¢1 + Amin (G2) lI tcc; 


t 
«V(x(t)) < (anes (H) + je?) Ix (|? + l e 1S) x1 (5)G1x(s)ds 


t-T 
t 
$ / e 9 33T(s)G,x (s) ds 
t-T 
1 : 
< (anes (H) + 5Bk cl?) [x (Ð|? + Amax (G1) xn, cs + Amax(G2)||Xller.c2- (2.2) 
These will be used in the paper. We also use the notation 
1 1 2 Amax(G1) Amax(G2) 
Ay := L— | Amax(H) + 2Bk|c|l^ | , Ago, = ———~., Ang, = ~~. 
17 icr (mall + ZB le?) «Anos = SEE Mao = eei 


The stability result formulated below depends, among others, on the positive definiteness of an auxiliary 
(3n + 1) x (3n + 1) matrix S. This matrix depends on matrices, vectors and constants used in (1.1)- (1.3), (2.1) 
and on a parameter v > 0, not yet involved. In other words, the elements of S are constructed using f, c, ¢2, 
k, v, b, c, A, B, D, H, G1, G2, formally 


S = S(B, c1, 62, K, V, b, c, A, B, D, H, Gi, G2). 
In its block form, the matrix S is defined as 
S44 S12 S15 S1 
Sh So S S 
Sis Sh S» S34 
Sh, Sh Si Sas 


with block-matrices 


S11 :- -ATH - HA - G4 - ATG2A, 

S1; := -HB - AT G;B, 

S13 :- -HD - AT G;D, 

Sie eue 7 (BAT + vr) c, 
S22 := eG, - B'G,B, 


$53 = =B'G2D, 


S24 := -B7 (62 + 56c) 
$33 = eG - p! G;D, 
Busse c + 3c) 


S44 = =b' Gab = Bc'b + n 


where I in 5,4 is the unit n x n matrix. 
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Theorem 1. If there exist non-negative constants D, v, positive constants ¢1, ¢2 and positive definite symmet- 
ric matrices H, G, and G2 such that the matrix S is positive definite, then the zero solution of system (1.1) is 
exponentially stable in the metric C° and the solution of the problem (1.1), (1.4) satisfies 


xi < (VA oO + Auc, Ilo + Vs Ils) e^. Q3) 


where t € [0, co) and fixed number y satisfies 


Amin (S) 
1 j* (2.4) 
Amax(H) + 5K |c| 


O0<y< y = min{¢,, 62, ô}, ô := 


Proof. The proof splits into two parts. In part i), an estimate of the total derivative of the functional (2.1) 
along solutions of the system (1.1) is computed. Part ii) deals with the exponential stability of solutions to 
system (1.1) in the metrix C°. 


i) Let us calculate the total derivative of the functional (2.1) along solutions of (1.1). Assuming t € [0, <œ) V 
(15,5 = 0, 1, ... we have 
Z V (x(t)) = [Ax(t) + Bx(t - T) + Dx(t - T) + bf (o ()| T Hx (t) 
+x" (t) H [Ax(t) + Bx(t - T) + Dx(t - T) + bf (o (t))] 
+ xT (t)G4x(t) - e *xT (t - T) Gix(t - 7) 
+37 (t)G5X(0 - e "x! (t - T) G3X(t- 7) 


t t 
-o1 f 699 (9) Grxls)ds - e; [ estt G2x(s)ds 


t-T tt 
+ Bf (o (£)) c? [Ax (t) + Bx (t - T) + Dx (t - T) + bf (o(0)]. 
Before writing this derivative in a form with the matrix S, we rewrite it using system (1.1) as follows 


£ V (x(t)) = [Ax(t) + Bx(t - 1) + Dx(t — 1) + bf (c (t))} T Hx (t) 


+x! (t) H [Ax(t) + Bx(t - T) + DX(t - 1) + bf (o (0)] 
«x! (0Gix(t) 
- e Ty? (t-t) Gax(t - 7) 
+ (Ax (t) + Bx (t - T) + DX (t — T) + bf (o ()). G2 
- (Ax (Ë) + Bx (t — T) + DX (t - 7) + bf (o (©) 
- ef?" (t - T) GaX (t - T) 


t t 
sea [ 699 endas - e; [ estt (9 G2x(s)ds 


t-T t-T 
+ Bf (o (£) c7 [Ax (t) + Bx (t - T) + DX (t - T) + bf (o(0)]. (2.5) 


Define an auxiliary vector 


X(t) := (rO x(t-1) Xx'(t-T) f(o()) . 
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Now, all the non-integral terms in the right-hand side of the last expression can be written using this vector 


and the matrix S. Omitting technicalities, it can be verified term by term that the identity 
[Ax(t) + Bx(t — T) + Di(t- T) + bf (o (£)]. Hx (t) 
+x! (t H [Ax(t) + Bx(t - T) + DX(t - T) + bf (o (0)] 
«x! (0Gix(t) 
- e Ty? (t — 7) Gax(t - 7) 
+ (Ax (t) + Bx (t - 7) + DX (t — T) + bf (a(t)? G2 
- (Ax (Ë) + Bx (t — T) + DX (t - 7) + bf (o (©) 
- e&t (t - T) GaX (t- T) 


+ Bf (o (©) c? [Ax (t) + Bx (t - T) + DX (t - 1) + bf (o (£))] 


- -X(0sx'(t) -v («co = if CW) f (c(t) 


holds. Then, (2.5) can be written as 


d v (x(0) =- XOXO - v G (0- if («») f (o (6) 


t t 


-a fox" (9) Gaxls)ds -ça f estt s) G2x(s)ds. 


t-T t-T 
Due to (1.3), 
-v (vo E («) f (o (5) « 0. 


As S is positive definite, we have 


FV (x(0) < -Amin (9) (KOJ? + [xt D|? + E-D +? (oO) 


t t 
-çı / e S15) xT (s) Gi x(s)ds - c; T e $357 (s) Go x(s)ds 
t-T t-T 


< - Amin (S) Ix (Ol? 


t t 


-çı / e S19) xT (s) Gy x(s)ds - c; I e S29) xT ts) Gy x(s)ds. 


t-T tt 


(2.6) 


(2.7) 


(2.8) 


ii) Prove the exponential stability of the solutions to system (1.1) in the metrix C°. Below we analyze two cases 


dealing with ç > 6 and ç < 6 where ç :- min{¢1, ¢2} > 0. 
ii;) Let ç > 6. Then, from inequality (2.2), we get 


t 


t 
- (Anas (H) + ipte) Ixco[ < -V(x(t)) + I e $C SyT(s)Gix(s)ds + Í e S20) 4T (s)Gzx (s) ds 


t-T t-T 


and 
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1 


- [x (9|? < ic 
Amax(H) + 3 Bk |c| 


V (x(t) 


t t 
" e 9U3)XT (s) Gi x(s)ds+ | e 9957 (s) G5x(s)ds | . 
1 2 
Amax(H) + > Bk [et X, 


cT t-T 


Inequality (2.8) yields 


1 
1 
Amax(H) + 5, Bk |c]? 


t t 
S 1 eS) xT (s) Gix(s)ds + | e € 9x (s) G5x(s)ds 
1 2 
Amax(H) + 3 Bk |c] 


A V (x(0) < Amin (S) V (x(b) 


t-T t-T 


t t 
-çı ] EEO Gix(s)as— c; Jamra G»X(s)ds 
t-T t-T 


or 


t t 
A v (x) < -6V(x(t)) - (61-8) [ TGs - (62-8) f e 9i! Gsxoas. 
t-T t-T 


We get, 
iy (x(t)) < -ôV (x(t)). 


Integrating this inequality over (0, t), which is a correct operation because the set of isolated points {t = 
TS, S =0,1,...}is countable, for t > 0, we obtain 


V (x (0) < VOe. (2.9) 


ii;) Let ç < 6. Rewrite the right-hand part of the inequality (2.2) as follows 


t t 
= ] $t eas = [ext @erxsds < -V(x(t) + (Anas (H) + jh?) Ixco |? . 
t-T t-T 


Now, using (2.8), we get 


F VECE) < -Amin (S) Ix (OP 
t t 
-ç | e9€*9XT (s) Gix(s)ds -ç | e "9x7 (s) Gaž(s)ds 
/ / 
< - Amin (S) Ix (D^ +ç (-v« (£) * (Ama +5k ie?) Ix oP) 
=- GV G0) - (Amin (S) - c (Ama) + ZBK eÈ) ) e 


=- çV (x(£)) - (6 - ç) (Ana + 5 Bk cl?) Ix (t)|? . 
Then, " 
ac V C) < -eV (x) 
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and integrating this inequality over (0, t), we obtain (as in ii4)) 

V (x(t) < V (p(0) e *', t2 0. (2.10) 
From (2.9) and (2.10) we deduce that, in both cases ii4), ii;) considered, we have 

V (x (0) < V (o(0)) e, t2 0. (2.11) 


Utilizing (2.2) and (2.11), we get (the below inequalities hold for t > O as well) 
Amin (H) [x(£)|^ + Amin (G1) llxll7 t,ç + Amin (62) Iži, < V (O) < V (9(0)) eo 
1 2 2 2 E) -yt 
< | { Amax (H) + 5 Bk |c] I9 (0)|^ + Amax (G1) lllo, c. + Amax (G) lllo, c.c | € 
hence, it follows that 
1 1 " z 
Ix(O|? < X5 [CE sie?) le 0) |^ + Amax(G1) lllo, c Amos (62) ls. e~, 


Now, after a simple estimation of the right-hand side, 


Ix (| < Folp, pe”? 


FolP, Q) :- (v Àn |9(0)| + \/Anc, IIXllo,r4¢ + \/ Ane, Ilo.) . 


Put no = 7/2. Inequality (2.3) is proved. The zero solution of (1.1) is C? exponentially stable within the mean- 
ing of Definition 2 since (1.6) holds. 


where 


3 Exponential C! Stability and Absolute Stability 


To prove the exponential stability in metric C1, we need the following lemma providing a formula transform- 
ing neutral system (1.1) to a delayed one of the Volterra type. 


Lemma 1. Let m € {1,2,...} be fixed. Then, for t € ((m - 1) v, mr), the solution of the problem (1.1), (1.4) 
satisfies the equation 


m-1 m-1 
x(t) = D™x(t - mr) + Ax(t) + 5 Dit (DA + B)x(t - ir) + D™1Bx(t - mr) + Y D'bf (o(t-ir)). (3.1) 
i-1 i-o 


Proof. We will prove the lemma by induction. Let x(t) be the solution of problem (1.1), (1.4). Then, for m = 1, 
the lemma holds because formula (3.1) coincides with initial system (1.1). Let us recall that, if we deal with a 
sum of the form X ..., then, by customary definitions, we set such sum equal to zero. This concerns the 
first sum in (3.1). Assuming that the conclusion of the lemma holds for m = j- 1 where j > 2 is a fixed natural 
number, we will show that, then, it holds for m = j as well. For m = j- 1, t € ((j-2)7, (j- 1)7), we get from (3.1) 


j-2 j-2 
X(t) = D'!x(t - G - 1)r) + Ax(t) + XO D™ (DA + B) x(t - it) + D'^Bx(t - (j- 1)r) + X` D'bf (o (t - it). 
i=1 i=0 


(3.2) 


If m =j and t € ((j - 1)r, jT), then, using (1.1), we get for the term x(t - (j - 1)7) in (3.2): 


x(t - (j - 1)r) = DX(t - jr) + Ax(t - (j - 1)t) + Bx(t - jt) + bf (o (t- G -1)r)). (3.3) 
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Finally, substituting (3.3) into the right-hand side of (3.2), we derive 


X(t) =D!" [Dx(t — jr) + Ax(t - (j - 1)t) + Bx(t - jr) + bf (o (t - G - 1)7))] + Ax(t) 


j-2 j-2 
+ XC D™ (DA + B)x(t - ir) + D^ Bx(t - (G - 1)t) + X` D'bf (o (t - ir) 

i-1 i-o 

n j-1 H 

-D'x(t - jr) + Ax(t) + X` D'* (DA + B)x(t - ir) 
i=1 
H jl H 
+ D'' Bx(t - jt) + S D'bf (o (t - it), 
i-0 


i.e., the formula (3.1) holds for m - j as well. 


Theorem 2. Let the hypotheses of Theorem 1, where y is fixed, hold and let 
0 < |De « 1. (3.4) 


Then, the zero solution of the system (1.1) is exponentially stable in the metric C1 and the solution of the prob- 
lem (1.1), (1.4) satisfies 


‘ B 8 " E 
HOME (Mvano DD Velo + Mu llello.s.c, * lllo + My/ Ane, lel... ) e, (3.5) 


where t € [0, œ) \ {Ts,s =0,1,...} and 


|DA + B| + |D| |b| |c| k 


M :5|A|* |b||c| k 
IAI + Ibl c| e. 


Proof. By Theorem 1, the zero solution of the system (1.1) is exponentially stable in the metric C°. We show 
that the zero solution of (1.1) is also C! exponentially stable. As it follows from formula (3.1) in Lemma 1, for 
t c ((m-1)t, mr), m=1,2,..., we have 


m-1 m-1 
Iž (Ð| < IDI" ||Xllo,e+ Al Ex ()|+|DA + B| Y ^ [D|* |x (t - iv) -LD]"" [B] lixclo +X 7 |DI' [b] If (o (t - in). 
i-1 i-0 


From (1.3) and (1.2) we derive |f (o (t - ir))| < k |o (t - ir)| < k|c||x(t - ir)| and, therefore, 
m-1 -— 
Ix (£) < IDI” [x]lo, + Al + |b| k lcl] Ix (| S 7 LDI [IDA + BI + [DI Ib|K]c]] |x (t — iv)] D] [BI lixllo,r - 
B (3.6) 
Let us estimate |x (t)| and |x (t — it)| in (3.6) by formula (2.3). We obtain 


(ÐI < IDI" lollo, + IDI IBI lollo, 
+ IIAL + bik ell (VA [O00] + s Ils. + es lose) ^ 
m-1 m 
+ [|DA + B| + |D] |b] k |c|]|D| ' - (S ren 
i-1 


(Vs let0)| + Au, Illo sss + / Ano, lel.) eth? 


- DI" lllo, + IDI” IBI lello, 
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+ [IA] + |bIk\cl] (Via |e(0)] + /Ane, lIPllo,r.¢, + Anc: Ilona) gat 


Á |DA + B| + |D| |b| k |c| [: ipi Quse] 
1 - |D| e? 


(Vu et] + ync: llo. s llos) 7. (3) 
Since t € ((m - 1)r, mr), then |D|" < exp((t/T)1n |D|). We have 
, z . B t 1 
DI" lli. + ID" Bl Plog < lei. + B Ilo, | exo (- tm i.) 


and, by (3.4), 


] dE" | 
IDI” llolio,s + IDI"? IBI rell. < [lllo + I5. lelo, | 77. 
ID] 


From (3.7), we get 


: : B _ 
HORE lilo. B lo eon 


" ja pices [DA + B| + |D| |b| k |c| n peer] 


1- |Diev/2 


(Vu eto s, lowe s llore). G8) 


Then, for any admissible m, inequality (3.8) with the term 


-|D te Date 


omitted, yields 
O| < Fil, pe”? 
where 
B 
31(9, p) = My/An |p(0)| + lllo, + My/Ane: Plor, + lllo, + M Ano, lllo, - 


Inequality (3.5) is proved. The zero solution of (1.1) is C! exponentially stable within the meaning of Defini- 
tion 3 where 54 = y/2 and (1.7) holds. 


The following theorem on absolute stability of system (1.1) within the meaning of Definition 4 directly follows 
from Theorem 1 and Theorem 2. 


Theorem 3. If the hypotheses of Theorems 1 and 2 hold, then the zero solution of the system (1.1) is absolutely 
stable. 


4 Example 


We will investigate system (1.1) where n = 2, T = 1, 


is -3 -2 , B- O 0.74035 D- 0.1 0 
1 0 0.74035 0 0 0.1 


T 
bev- (o. 01 -0.01) . o(t) = c' x(t) = 0.01xi(t) - 0.01x>(t), 
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i.e., the system 


X4(t) = 0.1% (t — 1) - 3x1(t) - 2x2(t) + 0.74035x2(t - 1) + 0.01f(o(t)), (4.1) 
Xx(t) = 0.1x2(t - 1) + 1xi(t) + 0.74035x4(t - 1) - 0.01/(o(t)), (4.2) 


with initial conditions (1.4). Set c1 = 0.01, c? = 0.001, k = 1, v = 0.1, f = 0.01 and 
0.5 0.1 0.06 0 0.67 0.43 
z = , H= ; 
e ( 0.1 02 | va ( 0 0.06 ( 0.43 0.67 
Part of the below numerical computations was performed by MATLAB. For the eigenvalues of matrices G4, 


G2 and H, we get Anis (G1) = 0.1697, Amax(Gi) = 0.5303, Amin(G2) = Amax(G2) = 0.06, Amin(H) = 0.24, 
Amax(H) = 1.1. The matrix S takes the form 


2.0600 1.5000 -0.3627 -0.3627 -0.0490 -0.0490 -0.0003 
1.5000 1.2800 -0.4960 -0.2295 -0.0310 -0.0670 0.0042 


-0.3627 -0.4960 0.4671 0.0990 O0 -0.0044 0.0005 
S= | -0.3627 -0.2295 0.0990 0.1651 -0.0044 O0 -0.0005 
-0.0490 -0.0310 O0 -0.0044 0.0593 O0 -0.0001 
-0.0490 -0.0670 -0.0044 0 Oo 0.0593 0.0001 


-0.0003 0.0042 0.0005 -0.0005 -0.0001 0.0001 0.0999 


and Amin(S) = 0.00000407. Because all eigenvalues are positive, matrix S is positive definite. Further we have 


EN 5Bk cl?) + 
_ Amax(61) + 
Aue, us) = 2.2095, 
_ Amax(G2) = 
AuG, "AED. v 0.25, 
5 -— Amin (S) - 3.7026 - 10%, 


Amax(H) + 5 Bk Ic? 
y -min (61, 62,0] = ô 
Moreover 
|A| 2 3.7025, |B| = 0.74035, |D| 0.1, |DA * B| 20.9186, |b| =|c| = 0.0141. 


Because 
y «-2r !In|D| = 4.6052 


we set y := 7" in (34). Then M = 4.7233. All hypotheses of Theorem 1 and of Theorem 2 are satisfied and, 
consequently, the zero solution of (4.1), (4.2) is exponentially stable in the metric C? and in the metric Ct. 
Finally, from (2.3), (2.4) and (3.5), it follows that the inequalities 


IKI < (Vs leC)| + Ane, llello,c.c, * Anc leo.) gu 

i (v4.5834 |9(0)| + V2.2095 || Pllo,1,0.01 + V0-25 | Pllo,1,0.001) g 90000018516. 
2 B 
soi < (t/a tO] + B los + MA llis 


E 5 —-—'t/2 
+ |dllos + My/Ane, Ilona ) et 


- (4.7233 4.5834 |P(0)| + 7.4 lllo, + 4.7233 V2.2095 |IPllo,1,0.01 


DE GRUYTER Josef Diblik et al., Absolute Stability of Neutral Systems with Lurie Type Nonlinearity —— 737 


+ |Illo,1 + 4.7233 0.25 loloa: ) B.V EN 


hold on [0, e») and on [0, e9)Ms, s =0,1,... ) respectively. By Theorem 3, the zero solution of the system (4.1), 
(4.2) is absolutely stable. 


Remark 1. System (4.1), (4.2) is considered in [5, 38, 44, 56] with omitted nonlinearity (f = 0), the same 
matrices A, D and with the matrix 
B-a O 1 , 
1 0 


where a is a constant. In [44] asymptotic stability is established for |a| < 0.4 while, in [38], for |a| < 0.533 
and, in [56] as the authors state, for |a| « 1. In these papers no estimates are given related to the convergence 
of solutions or their derivatives. Exponential stability with explicit formulas estimating the convergence of 
solutions and their derivatives is given in [5] for |a| < 0.6213. The present example gives exponential stability 
and explicit formulas estimating the convergence of solutions and their derivatives for |a| < 0.74035. 


5 Comments and Conclusions 


The direct Lyapunov method is widely used for the stability analysis of functional differential systems. Its suc- 
cessful application substantially depends on an ingenious construction of Lyapunov-Krasovskii functional, 
suitable for the problem considered. As a rule, in the investigation of nonlinear neutral type control systems, 
functionals are used in the form of the sum of a quadratic form of terms on the right-hand side of the given 
system, a quadratic form of a *prehistory" and an integral of the nonlinearity. For example, in [27, 51], in the 
case of system (1.1), the following one is used: 


t o(t) 


VG) = [x(t) - Dx(t - 0] H [x(t) - Dx(t - 0] + / eS) XT (s)Gx(s)ds + B f f(s\ds (5.1) 
0 


t-T 


where H and G are positive definite matrices, c > 0 and B > 0. However, applying functional (5.1) or any such 
functional type, does not result in an estimate of the norm |x(t)| of solutions, but rather in the integral norm 
IX \lt,7,¢ Many authors use method of functionals or other methods to state the fact of asymptotic stability only 
because their approach does not give an explicit estimate of the norm |x(t)| or |x(t)| (such as formulas (2.3), 
(3.5) in Theorem 1 and Theorem 2). We refer, e.g., to [7, 18, 23, 37, 56]. Functional (2.1) is more general than some 
of the previously ones used to investigate neutral equations. For example, functionals used in [5, 35, 52, 53] 
are, if the systems considered are reduced to the form (1.1), its particular cases. 

Let us discuss some specific properties of our approach. By an appropriate choice of the coefficient v in 
the blocks S43 and S44 of the matrix S, in some cases, the positive definiteness of S can be achieved. Such a 
parameter is involved because of the sectorial condition (1.3) used in the proof (inequality (2.7)). The block- 
matrix S4; contains the expression -ATH - HA. By the well-known result on Lyapunov matrix equations, if 
system x(t) = Ax(t) is asymptotically stable, then, for any positive definite symmetric matrix Q, there exists 
a unique positive definite symmetric matrix H such that -ATH - HA - -Q. This fact can be, in principle, 
used when the derivative of the functional V along trajectories is estimated (see (2.5), (2.6)). We did not use 
such a possibility preferring consideration with the matrix S giving more general results. The parameters c1, 
¢2 in the exponential multipliers exp(-c1(t - s)), exp(—¢2(t — s)) in the definition of V(x(t)) by (2.1) have an 
impact on diagonal block-matrices S22, S33 and, moreover, as suggested by formulas (2.4), (3.5), on the speed 
of convergence of the vanishing solutions. As it follows from the proof of Theorem 1, part ii2), formula (2.10) 
(and computations above this formula), the assumption c4 > 0, c? > O cannot be replaced by çı > 0, ¢2 > 0. 
If c4 = Oor if c? = O then, proceeding by the scheme of the proof, we are not able to prove an exponential 
stability and construct explicit formulas estimating solutions and their derivatives. 
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In Theorem 2 we assumed |D| # 0. This is a quite natural assumption because if |D| = 0, then, by (1.5), 
all eigenvalues of the matrix DTD are zeroes. Consequently, D = 0, where @ is n x n zero matrix and (1.1) 
becomes a delayed rather than neutral system and can be investigated with a functional defined by (2.1) if 
the integral with kernel defined by the derivatives is omitted. In many investigations when neutral systems 
such as (1.1) are considered, assumption |D| « 1 is necessary. In the present paper, this inequality follows 
from formula (3.4) in Theorem 2. In Theorem 1, a block-structured matrix S is applied. Although the latter 
assumption is not explicitly formulated (and is not explicitly used in its proof), it is probably implied by the 
positivity of S. 

It is an open question, whether a functional Vp(x(t)) = e"*V(x(t)), where p is a suitable real parameter 
and V(x(t)) is defined by formula (2.1), can be used to estimate solutions to system (1.1) in non-stable cases. 
Further progress can be achieved as well if the approach used is extended to systems (1.1) with the sum of 
nonlinearities X bifi(o;(t)), o;(t) = e x(t) of a Lurie type instead of only one such nonlinearity. Another 
challenge is a possible generalization of the results to systems with variable matrices and delays and to 
fractional differential systems as well. 
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